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ABSTRACT 

The probability distribution function of column density (N-PDF) serves as a powerful tool to characterise the various physical pro¬ 
cesses that influence the structure of molecular clouds. Studies that use extinction maps or H 2 column-density maps (N) that are 
derived from dust show that star-forming clouds can best be characterised by lognormal PDFs for the lower N range and a power-law 
tail for higher N, which is commonly attributed to turbulence and self-gravity and/or pressure, respectively. While PDFs from dust 
cover a large dynamic range (typically N cm“^ or Av~0.1-1000), PDFs obtained from molecular lines - converted into H 2 

column density - potentially trace more selectively different regimes of (column) densities and temperatures. They also enable us to 
distinguish different clouds along the line of sight through using the velocity information. We report here on PDFs that were obtained 
from observations of ^^CO, ^^CO, C^^O, CS, and N 2 H^ in the Cygnus X North region, and make a comparison to a PDF that was 
derived from dust observations with the Herschel satellite. The PDF of is lognormal for Av~l-30, but is cut for higher Ay 
because of optical depth effects. The PDFs of and are mostly lognormal up to Av~l-15, followed by excess up to Av~40. 
Above that value, all CO PDFs drop, which is most likely due to depletion. The high density tracers CS and N 2 H^ exhibit only a 
power law distribution between Av~15 and 400, respectively. The PDF from dust is lognormal for Av~3-15 and has a power-law 
tail up to Av~500. Absolute values for the molecular line column densities are, however, rather uncertain because of abundance and 
excitation temperature variations. If we take the dust PDF at face value, we ‘calibrate’ the molecular line PDF of CS to that of the 
dust and determine an abundance [CS]/[H 2 ] of 10"^. The slopes of the power-law tails of the CS, N 2 H+, and dust PDFs are -1.6, -1.4, 
and -2.3, respectively, and are thus consistent with free-fall collapse of filaments and clumps. A quasi static configuration of filaments 
and clumps can also possibly account for the observed N-pdfs, providing they have a sufficiently condensed density structure and 
external ram pressure by gas accretion is provided. The somehow flatter slopes of N 2 H^ and CS can reflect an abundance change 
and/or subthermal excitation at low column densities. 

Key words. ISM: clouds, ISM: extinction, ISM: abundances, ISM: molecules, ISM: structure 


1. Introduction 

Probability distribution functions (PDFs) form the basis of many 
modern theories of star formation (e.g. Krumholz & McKee 
I2005I Hennebelle & Chabrier [2008l Federrath et al. l2008ll2010i 
Padoan & Nordlund I201H Federrath & Klessen I2Q121120131 
Hopkins 12013b . and are frequently used to characterise prop¬ 
erties of the interstellar medium in simulations (Klessen 120001 
Vazquez-Semadeni & Garcia [200T1 Burkhart et al. 120131 Ward 
et al. 12014b . In summary, a PDF is defined as the probability 
of finding gas within a column-densit}Q range [N, N-\-dN]. We 
define r] = ln(N/(N)), and the quantity Pr^ip) then corresponds 


Correspondence to: nschneid@phl.uni-koeln.de 
^ The column density N can be expressed as visual extinction Ay 
with N(H 2 ) = Ay 0.94x10^^ cm“^ mag“^ (Bohlin et al. ll9781) . 


to the PDF of p. For more details on the basic definitions, see 
Schneider et al. (I2015ab . 

Observationally, PDFs were first obtained from near-IR ex¬ 
tinction maps (Lombardi et al. 120081 Kainulainen et al. 1 20091 
Froebrich & Rowles [2010b . With the advent of Hersche^ it is 
now possible to determine dust column-density maps that cover 
a very large dynamic range from Ay<1 up to a few hundred Ay. 
Studies using this type of map reveal large variations in dust 
PDFs. Earlier findings (Kainulainen et al. 2009) of a simple 
lognormal plus power-law tail distribution for star-for ming re¬ 
gions are co nfirme d (e.g. Schneider et al. 12013112015a[ Stutz & 
Kainulainen 12015b . but more complex shapes are also detected 
(Hill et al. 120111 Schneider et al. 12012112015cl Russeil et al. 

^ Herschel is an ESA space observatory with science instruments pro¬ 
vided by European-led Principal Investigator consortia and with impor¬ 
tant participation from NASA. 
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Fig. 2. Line integrated (v=-10 to 20 km s maps (in K kms of molecular emission in Cygnus X North. The lower level is zero 
for all maps. The blue crosses in the N 2 H^ map indicate submm-continuum sources (Motte et al. 120071) . 
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(either as a phase transition between clumps and interclump gas 
or stellar feedback). 

Only a few studies attempt to make PDFs from molecular 
line observations. A/^(H 2 )-PDFs of or (Goldsmith et 
al. 120081 Wong et al. 120081 Goodman et al. 120091 Lo et al. 120091 
Carlhoff et al. 120131 Schneider et al. 12015bl) are often clipped 
at a certain Ay threshold because the lines can become optically 
thick and CO is depleted. Taking W43 as an example, this thresh¬ 
old for ^^CO 2^1 lies between Ay = 30 (for an assumed T^y^lO 
K) and -100 (for Tex=12 K and correction for high opacity). 
Assembling a PDF from one or several molecular line tracers 
with high critical densities such as CS, HCO^, HCN, N 2 H^ etc. 
has not yet been tried owing to the lack of extended maps. In 
this study, we use our large dataset of Cygnus X in the molecular 
lines of i^CO, ^^CO, and C^^O 1^0, CS 2^1, and N 2 H+ 1^0 
that were obtained with the Five Colleges Radio Astronomy 
Observatory (FCRAO), and dust column-density maps of the 
Cygnus X North region (Hennemann et al. 120121 Schneider et 
al. I2015dl) to produce PDFs from H 2 column density that were 
obtained from molecules and dust. We discuss their properties 
but also point out the large uncertainties and observational diffi¬ 
culties related to molecular line PDFs. 


2. Observations 


Fig. 1. 1^0 contours of line integrated (v=-10 to 20 km 

s“^) emission with levels 63.5 to 381 K kms“^ in steps of 63.5 
K kms“^ overlaid on a Herschel dust column-density map at 36" 
resolution (that corresponds to 0.41 pc at a distance of 1.4 kpc). 
This map was corrected for an average foreground contamina¬ 
tion of Ay= 5. Prominent features, such as W75N or the DR21 
ridge, are indicated in the plot. DR17, 22, and 23 are H ii regions. 


120131 Alves de Oliveira [20141 Tremblin et al. 120141) . The out¬ 
come of these works is that (i) all clouds have a distribution for 
low Ay that is best described by a lognormal, (ii) a second peak 
at low Ay can emerge in the case of external pressure, which rep¬ 
resents the compressed shell, and (iii) a single or double power- 
law tail is found for star-forming regions, and it is proposed that 
this is related to gravitational collapse and/or external pressure 


Molecular line data 

We use data from a molecular line survey of (115.3 GHz), 
(110.2 GHz), 1^0 (109.8 GHz), CS 2^1 (98.0 
GHz), and N 2 H^ 1^0 (93.2 GHz), which were collected with 
the FCRAO 14m radiotelescope between 2003 and 2006. The 
data were obtained using the 32 pixel Second Quabbin Optical 
Imaging Array (SEQUOIA) in an on-the-fly (OTF) observing 
mode. The beamwidth of the FCRAO at 93 GHz is -48" and at 
110 GHz -45", while the main beam efficiency at that time was 
0.48. For the present paper, we use maps on a 20" grid. The data 
have a mean 1 cTrms I'ms noise level of -0.4 K per 0.06 km s“^ 
channel on a Tmb temperature scale that we use in this paper. In 
the C^^O map, the noise level and a slightly stripy appearance, 
resulting from the on-the-fly mapping mode, becomes apparent. 
However, this does not affect the PDFs because they are assem¬ 
bled from pixels above the ScTrms level. For further details, see 
Schneider et al. (IMTOllMTTI) . 
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Column-density maps from Herschel 

For this study, we employ a column-density map at 36" angular 
resolution (regridded to 20" to match the fully-sampled FCRAO 
maps) that were obtained as part of the HOBYS {Herschel imag¬ 
ing survey of OB Young Stellar objects) key program (Motte et 
al. 120101) . This map was obtained by an SED (Spectral Energy 
Distribution) fit to the PACS (Poglitsch et al. l20101) 160 yum, and 
the SPIRE (Griffin et al. 120^ 250, 350, and 500 jim wave¬ 
lengths observations. More details can be found in Appendix A 
and a full map of the column de nsity of the Cygnus X North re¬ 
gion is shown in Schneider et al. (I2015db and Bontemps et al. (in 
prep.). A smaller cutout of the DR21 ridge at 25" resolution is 
presented in Hennemann et al. (120121) . 

3. The structure of the Cygnus X North region 

Cygnus X is one of the richest star-formation regions in the 
Galaxy, containing the prominent OB-association Cyg OB2, 
with ~50 0-stars (see Reipurth & Schneider 120081 for a review). 
Erom large-scale ^^CO 2^1 (Schneider et al. 120061) and ^^CO 
1^0 (Schneider et al. 120071) surveys, we derive a mass of a few 
10^ M© for the whole molecular cloud complex that is divided 
into the Cygnus X North and South regions. These studies show 
that the majority of the molecular clouds in the complex are lo¬ 
cated at a common distance of about 1.4-1.7 kpc. A distance 
of 1.4 kpc, which is derived from maser parallax (Rygl et al. 
|20I3, is now commonly accepted. The clouds in the Cygnus X 
North region are actively forming stars and contain more than 
100 massive pre- and protostellar dense cores (Motte et al. 120071 
Bontemps l20101 Csengeri et al. 1201 II Hennemann et al. 2012). 

Eigure shows the dust column-density map of Cygnus X 
North with contours of 1^0 overlaid onto it. The 
emission is integrated over the velocity range -10 to 20 km s“\ 
which confines all star-forming molecular clouds that are di¬ 
rectly associated with the Cyg OB2 cluster. Outside this veloc¬ 
ity range, there are no other clouds along the line of sight, al¬ 
though in other regions of Cygnus X, clouds from the Perseus 
arm appear at velocities around -40 km s“^. As demonstrated 
in Schneider et al. (120061 see their Pig. 3), there are mainly 
two velocity coherent cloud complexes in Cygnus X North: the 
DR21 ridge-DR22-DR23 clouds at v=-10 to 1 km s“^ and the 
W75N-DR17 clouds at v=7 to 20 km s“^ However, and 
dust also trace more diffuse emission, mainly arising from the 
‘Great Cygnus Rift’, a nearby (0.6-1 kpc) region of obscuration, 
which is identified in optical images, and which is not associ¬ 
ated with Cyg OB2. As shown in Schneider et al. (12007b . the 
extinction due to the Rift is of the order of Ay-5-10 and mainly 
arises between v=6 to 20 km s“^ Other authors determine val¬ 
ues between Av=2-5 (Sale et al. 12009b . Av=2.5-7 (Drew et 
al. f2008b . Av=5. 5-7.5 (Wright et al. l20TQb . and Av=2.6-5.6 
(Guarcello et al. 12012b . An extinction around five is also the 
lowest emission level in the Herschel dust map that is analysed 
in this paper. The Rift emission is barely detected in and 
not in C^^O, CS, and N 2 H^. This means that we expect only 
the low column-density range of the PDP to be affected. 
The higher column-density power-law tail is, in any case, not 
affected. 

Pigurej^shows the same cutout of Cygnus X North as Pig.[2 
in four different line tracers. The optical depth of the lines de¬ 
creases from left to right and the critical density increases. While 
the CO 1^0 lines have a low critical density ricr of ~10^ cm"^, 
the CS 2^1 and the N 2 H^ 1^0 lines have ricr of 1.3x10^ cm“^ 
and 6.lx 10"^ cm"^, respectively, for a temperature of 10 K and 
thermal excitation (see. e.g. Plower 119991 Daniel et al. 120051 
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rj = LN(N/<N>) 

Fig. 3. PDP derived from the foreground-corrected Herschel 
column-density map. The left y-axis gives the normalized prob¬ 
ability density the right y-axis the number of pixels per log 
bin. The upper x-axis is the visual extinction and the lower x- 
axis the natural logarithm of the normalized column density. The 
green curve indicates the best lognormal fit to the low column- 
density distribution, the red line displays a linear regression 
power-law fit to the high column-density tail. The dispersion of 
the fitted PDP (cr^), the slope 5- and its error, and the exponent 
a of an equivalent spherical density distribution (C) and a cylin¬ 
drical density profile, such as for filaments (P), are indicated. 



Av 


Fig. 4. The original PDF of a model cloud (see text) with a log¬ 
normal core and a power-law tail is shown with a green dashed 
line. This cloud is ‘contaminated’ with a model cloud that has 
a lognormal core. The resulting PDF is shown as a blue dashed 
line. The red line then displays the ‘corrected’ PDF if a constant 
value of Ay=5 is removed. 


Shirley et al. 12015b . In addition, CO depletes in cold and dense 
gas while CS, and in particular N 2 H^, remain stable in this gas 
phase. Accordingly, the maps show selectively different density 
regimes of the gas. While 1^0 is still sensitive to lower- 
density gas and resembles the map, N 2 H^ is confined to 
the densest clumps. Basically all clumps visible in the N 2 H^ 
map contain cold, dense cores (starless and protostellar, indi¬ 
cated by crosses in Fig.[^ observed in mm-dust continuum with 
the MAMBO bolometer (Motte et al. 12007b . From the submm- 
dust observations, a typical size of 0.1 pc (0.7 pc) and density 
of 10^ cm“^ (10"^ cm“^) for the cores (clumps) was derived. This 
value is confirmed by a decomposition of our N 2 H^ map using 
the GAUSSCLUMPS algorithm (Stutzki & Glisten [T990b in the 
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velocity range -10 to 20 km s“\ which shows an average clump 
size of -0.6 pc. We only consider clumps larger than 1.5 times 
the beamsize (45"). 



Fig. 5. Probability distribution functions of H 2 , obtained from 
dust and ^^CO, ^^CO, and 1^0. All PDFs were con¬ 

structed using pixels above the 3cr noise level. 



Fig. 6. Like Fig. 4 but showing the CS 2^1 and N 2 H^ 1^0 
PDFs (converted into H 2 column density) in addition to the dust 
PDF. 


4. Probability distribution functions from dust and 
molecules 

4.1. The dust PDF 

As discussed in the previous section, the dust map may suf¬ 
fer from foreground contamination that mainly arises from 
the Cygnus Rift. Contamination always modifies the measured 
cloud PDF. Relative to the underlying PDF, the measured PDF 
becomes narrower, the peak shifts towards higher Ay, and the 
slope of the power-law tail steepens (Schneider et al. 12015^ . 
Following the method outlined in Schneider et al. (I2015al) . we 


can correct for foreground contamination by removing the con¬ 
stant Ay=5 value for the Cygnus Rift from the original column- 
density map. 

For a more realistic approach, we have to take into acount, 
however, that the foreground contamination will typically also 
show a distribution of column densities, usually also following 
a log-normal distribution. In a systematic study, Ossenkopf et 
al. (in prep.) show that the correction by a constant foreground 
subtraction still works quite well if the most probable contami¬ 
nating column density is used as the offset. This is i llustrate d in 
Fig.0 Following the description in Schneider et al. (I2015al) . we 
contaminate a model cloud with the parameters of the Cygnus 
PDF with a foreground cloud that is generated from a fractional 
Brownian motion (fBm) map. This has a power spectral index 
of 2.8, which matches the typically-observed spatial scaling be¬ 
haviour, and has a lognormal column-density PDF with a width 
(Trj = 0.45 and a contamination Ay= 5.0. This cloud is thus a 
more realistic representation of a foreground or background con¬ 
taminating cloud. Figure shows that the contamination pro¬ 
duces a distortion in the PDF but the simple correction obtained 
by subtracting a constant Ay=5.0 value fully recovers the PDF 
tail and recovers the PDF peak position approximately. 

Even for such a relatively strong contamination, we can thus 
recover the main PDF properties by subtracting a constant off¬ 
set for the contamination. Figure shows the N-PDF obtained 
from Herschel dust continuum data, corrected for a contamina¬ 
tion of Ay=5. This dust PDF shows the typical shape of PDFs 
obtained for star-forming regions: a lognormal distribution for 
low column densities followed by a power-law distribution be¬ 
tween Ay-10-15 and a few hundred Ay. Above Ay-100, some 
excess in the PDF is observed that could constitute a second 
power law, when compared to a simple power law. Such an ex- 
cess/fiatter power-law tail was recently found in some massive 
GMCs (Schneider et al. l2015cl) and was interpreted as a possible 
result of internal stellar feedback. The same argument may hold 
here because only pixels in the DR21 ridge, W75N, and some 
small, isolated clumps - all associated with young stellar objects 
(YSOs) - contribute to this high column-density part of the PDF 
(see Fig.[^. 



Fig. 7. N-Probability distribution functions obtained from CS us¬ 
ing (i) different excitation temperatures (5, 10, 15 K) and an 
abundance of 4x10“^^, and (ii) at a temperature of 10 K but 
abundance values of 10“^^ and 10“^. The dust PDF is shown 
for comparison. 
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4.2. The PDFs from CO 

The PDFs of H 2 column densities N^o/(H 2 ), which were ob¬ 
tained from molecular lines, are shown in Fig. andtogether 
with the same dust PDF from Fig. Appendix B explains in 
detail how N^o/(H 2 ) was derived. In summary, we assume a 
constant excitation temperature of 10 K, LTE conditions, and a 
beam/source-filling factor of unity when determining all column 
densities. We note, however, that PDFs are always resolution- 
limited. The PDFs were constructed from pixels above the 3cr 
noise level and are normalized to the number of pixels from the 
dust map (ffpix^ust), i-e. we scaled the molecular line PDFs by 
the factor #pix^o//#pix^^^^. In Appendix C, we show the individ¬ 
ual PDFs without this scaling. Above the noise level of Av~l-2, 
the N„o/(H 2 )-PDFs of ‘^CO, ‘^CO, and C‘^0 can be fitted by 
a broad lognormal distribution with a peak at Ay-3-5, and by 
a power-law tail above Av=15 for C^^O. Also, the PDF 
shows some excess above Ay-15 (see Appendix C for a more 
detailed display). 

Considering all uncertainties in the determination of H 2 col¬ 
umn density from dust and CO, the lognormal part in the CO- 
PDFs corresponds quite well to the lognormal part of the dust 
PDF, which indicates that the low-density molecular gas is prob¬ 
ably well mixed with dust. The ^^CO PDF is a special case 
because it is contaminated by line-of-sight confusion from the 
Cygnus Rift at a distance of 600 pc, and thus shifts towards lower 
Ay to some degree. For Ay larger than ~15 (we note that the col¬ 
umn densities are correct only within a factor ~2), however the 
^^CO PDFs depart from the lognormal shape and the distribu¬ 
tion becomes ‘bumpy’ and falls off. This is the regime where the 
line no longer traces the dense gas because it becomes optically 
thick. In comparison to hydrodynamic simulations with radia¬ 
tive transfer (Shetty et al. 1201 IK our observed ^^CO PDF covers 
a lower column-density range. Their ‘Milky Way cloud’ (case 
(a) in Fig. 2) PDF has the same width but pea ks around 
Ay=10-15. More recent results of Sziics et al. (120151) . however, 
show PDFs that cover a very similar column-density range 
to the one we observed. 

The and PDFs extend to higher Ay (up to ~40). 
The line is optically thin (r ~0.1-0.3, see Fig. B.3) in 
most of the clouds and becomes marginally more optically thick 
(r ~1) only in the DR21 ridge. The line is optically thin 
everywhere with maximum values of r ~0.4 in the DR21 ridge. 
Optical depth effects, and not depletion, are the main reason for 
the different PDF shapes of ^^CO, ^^CO, and for higher 
column densities. The power-law behaviour for Ay>15 that we 
hnd for indicates that the whole column of gas is traced 
up to Ay~40. Because becomes marginally optically thin, 
we do not observe a clear power-law tail but only some ex¬ 
cess because the highest column densities are no longer traced. 

becomes even more easily optically thick so that the high 
column-density pixels are not traced at ah and the PDF appears 
lognormal. Depletion of ah CO isotopologues probably sets in 
for Ay larger than around 40 where ah PDFs fall off. The re¬ 
moval of CO from the gas phase typically happens at d ensitie s 
above 10"^ cm“^ at temperatures ~15 K (e.g. Blake et al. 119951) . 
Our value of Ay~40 (A= 3.8x10^^ cm“^) thus corresponds to a 
typical clump size d of -1.25 pc (d = A/(10'^cm“^)) in which 
CO depletes. These clumps are the ones that are well traced by 
N 2 H^ emission and the temperature map (Fig.[^ that we discuss 
in the next section. 

At the low Ay end, approximately in the range Ay=1-2, we 
observe some excess compared to a purely lognormal PDF dis¬ 
tribution for ^^CO and C^^O (Fig. C.l). The chemistry in this 


column-density regime is governed by (F)UV radiation via pho¬ 
todissociation (Ay<5) and fractionation (Ay<3) processes (e.g. 
van Dishoeck & Black [T988]) . The excess indicates that self¬ 
shielding of ^^CO and C^^O against UV-radiation is rather ef¬ 
fective. Otherwise, the abundances would be lower and the PDF 
would drop more signihcantly. On the other hand, in warmer re¬ 
gions (internally or externally heated) in the medium Ay range, 
CO desorption leads to a higher abundance of CO. These spatial 
abundance variations, which were observed for ^^CO in Perseus 
(Pineda et al. 12008b and Orion (Ripple et al. l2013b . have an obvi¬ 
ous influence on the N-PDF. However, they are difficult to assess 
(see next Section). 

4.3. The PDFs from CS and N 2 H^ 

The N^o/(H 2 )-PDFs of CS and N 2 H^ mostly cover the high 
column-density range and can be best htted by a pure power 
law (see Appendix C). However, ah absolute values need to be 
taken with caution because they strongly depend on the abun¬ 
dance values that were used for converting column densities into 
H 2 and on the excitation temperature. 

For N 2 H^, we adopt [N 2 H^]/[H 2 ] = 5x10“^^, a value deter¬ 
mined for massiv^ cloud cores (Pigorov et al. 12003b . and ap¬ 
plied it to one of the Cygnus X North cores (N63, Motte et al. 
2007) in the following way (Fechtenbaum et al., in prep.): Fevel 
population diagrams that use lines of the optically thin isotopo¬ 
logues 15 N 2 H+ and N^^NH^ were plotted to estimate their col¬ 
umn density. An isotopic ratio of 448 (Wilson & Rood 

11994b was used to derive the column density of N 2 H^. The mass 
of 44 Mo in a 30" beam (Duarte-Cabral et al. l2013b of the source 
enabled us to calculate the column density of H 2 and thus the 
abundance of N 2 H^. The resulting PDF covers a column-density 
regime of Ay-15 to 400, which corresponds very well to the one 
derived from dust (see Fig.[^. 

For CS, there is a large scatter in observed values for high- 
mass star-forming regions with [CS]/[H 2 ] = 10~ ^ to a few 
ti mes 10 ~^^ (Gianetti et al. 120121 Fi et al. 120151 Neufeld et 
al. 12015b . Fow-mass c ores h ave typical values of a few times 
10"^ (e.g. Tafaha et al. 12006b . Here, we use a value of 4x10“^^ 
(Fechtenbaum et al., in prep.) that was obtained using the op¬ 
tically thin isotopologue C^^S and applying the same method 
as outlined above. The PDF obtained with this abundance value 
covers the Ay range -20 to 350 and thus has a similar column 
density range as the PDF determined from dust, but shifted ver¬ 
tically to the latter. 

To investigate in more detail the effects of different excita¬ 
tion temperatures and abundances, in Fig. [7] we plot a whole set 
of PDFs for CS. We hrst vary the excitation temperature, using 
5 and 15 K as two extreme cases, and 10 K (average value de¬ 
termined from ah molecular lines and dust). The column density 
variation is signihcant, increases its value by more than a factor 
two between Tex=15 K and 5 K. Different abundance values ob¬ 
viously have a signihcant impact on the PDF as well. A lower 
(but unlikely) [CS]/[H 2 ] ratio of 10“^^ shifts the PDF to higher 
column densities with unrealistic values of a few hundred Ay . 
A higher abundance shifts the PDF towards lower column den¬ 
sities, and a value of 10“^ leads to a PDF that corresponds to 
the one obtained from dust. In this way, it is possible to ‘cal¬ 
ibrate’ the [CS]/[H 2 ] abundance using the dust PDF. However, 
this method should be treated with caution because the molecu- 


^ Lower abundances around 1-3x10“^^ were found in low-mass 
dense cores and da rk clouds (e.g. Caselli et al. 120021 Bergin et al. 1200 11 
Tafalla et al. 120061) . 
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Hg column density from Herschel 



Temperature from Herschel, Av=15 contour from NoH^ 



RA(2000) 


Fig. 8. Left: Foreground-corrected dust column-density map (as shown in Fig. 1) in greyscale. The black contour indicates the Av= 
15 level. The red contours show the same Av=15 level from the H 2 column-density map derived from N 2 H^, using an excitation 
temperature of 10 K and an abundance of 2x10“^. Right: Dust temperature map of Cygnus X North from Herschel with the same 
red Av=15 contour from N 2 H^. 


lar line PDF also depends on the excitation temperature, and the 
dust PDF has its own uncertainties (variable temperature along 
the line of sight, unknown dust opacity etc.). 


4.4. The high column density range in the maps 


The PDFs obtained from dust emission and N 2 H^ (and CS) both 
show a power-law tail starting around Av~15. Figure (left 
panel) illustrates the high spatial correspondence between the 
pixels in the dust and N 2 H^ maps that cover the same column- 
density range. The black and red contours indicate the Ay = 15 
level for the H 2 column density derived from dust and N 2 H^, 
respectively. On average, both contours cover the same areas, 
implying that the dust is well mixed with the dense and cold 
gas that is traced by N 2 H^. The geometry of the emitting re¬ 
gions is a mixture between spherical and filamentary. Even at 
very high column densities (indicated by the arbitrary chosen 
Ay=100 blue contour from N 2 H^), not only spherical clumps 
are outlined but also the elongated DR21 ridge. An overlay be¬ 
tween the dust temperature obtained with Herschel and the same 
Ay=15 contour from N 2 H^ (Fig.[^ right panel) shows that it is 
indeed mostly cold gas traced by N 2 H^ and dust. The clumps 
that are indicated by the red contour typically have a tempera¬ 
ture of 10-15 K, except those with internal sources (the DR21 
region in the DR21 ridge and DR23 at the southern border of 
the map). In these cases, the column density is underestimated 
so that the true power-law tails of both, the dust, and N 2 H^ PDF, 
which is shown in Fig.|^ are slightly fiatter. 


Species 



ac 

Aac 

a/ 

Aaf 

Dust 

-2.29 

0.03 

1.87 

0.03 

1.44 

0.03 

C^^O 1^0 

-2.57 

0.12 

1.78 

0.04 

1.39 

0.04 

CS 2^1 

-1.56 

0.02 

2.28 

0.02 

1.64 

0.02 

N 2 H+ 1^0 

-1.41 

0.03 

2.42 

0.03 

1.71 

0.03 


Table 1. Results from the power-law fit to the PDFs. 5’ and A a are 
the slope and its error, ac and af are the exponents of a density 
distribution with p oc r~^ with spherical and filamentary geome¬ 
try. 


5. Discussion and conciusions 

In the following section, we base our discussion on the column 
density maps and PDFs of CS and N 2 H^ derived using an exci¬ 
tation temperature of 10 K and abundance values of [CS]/[H 2 ] 
= 10“^ (to be consistent with the dust PDF) and [N 2 H^]/[H 2 ] = 
5x10“^^, respectively. 

5.1. Thermal and subthermal excitation regimes 

For a homogeneous medium with beam-filling of unity, the H 2 
column-density threshold N for thermally excited N 2 H^ and CS 
calculates from the clump size (d) and critical densities with 
N = d X Her- For d we a dopt 0 .7 pc, taken from dust submm 
observations (Motte et al. I2007I) . This value is close to the av¬ 
erage clump size of 0.6 pc we derive from our N 2 H^ clump 
decompositioij^ We thus obtain Ay=285 and 130 for CS and 

^ The submm continuum observations have an angular resolution of 
11", compared to the 45" for N 2 H^, and trace a much larger column 
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N 2 H^, respectively, assuming critical densities of 1.3x10^ cm“^ 
(6.1x10"^ cm“^) for CS (N 2 H^). Above these values, the emis¬ 
sion lines should be thermalized. However, the dumpiness of 
the gas leads to a lower threshold of thermalization. With an 
average density cm“^ of the clumps,which was deter¬ 

mined from submm continuum observations (Motte et al. 2007), 
we obtain a beam filling ((^)/^cr) of -10% for CS and -15% 
for N 2 H^. Both lines are thus already thermalized at lower Ay . 
For N 2 H^ we independently derive its thermal excitation by the 
hyperfine-structure line ratios. The line is mostly thermally ex¬ 
cited because we obtain excitation temperatures above 5 K for 
the emitting gas (Appendix C), which is above the typical value 
for subthermal excitation. 

CS has a higher critical density but depletes at densities 
above -10^ cm“^, which is consistent with our PDF that is cut off 
at Ay-100 (for the abundance [CS]/[H 2 ] = 10“^). On the other 
hand, its chemistry is less density dependent. Hence, we expect 
that CS is not thermalized at the lowest end of the distribution. 

5.2. Slopes of the N-PDFs 

Independent of the exact column-density regime that is covered 
by the PDFs from dust and CS/N 2 H^, they all show power-law 
tails with slopes 5' between -1.4 and -2.3 (see Table 1). The 
slope 5' of the power law can be converted into the exponent 
a of an equivalent density distribution p oc r~^ (Appendix D). 
For spherical geometry, representing a single core or a core en¬ 
semble (Girichidis et al. 12014b . we use ac=l-2ls (Federrath & 
Klessen l2Q13b . For singular polytropic cylinders, which portray 
filaments, we use af=l-l/s (see Toci & Galli (12015b for regu- 
lar pol ytropic cylinder models and Fischera (12014b and Myers 
(12015b for the corresponding N-PDFs). The values we obtain for 
a (ac or af) vary between 1.4 and 2.4. In the simplified picture 
of the free-fall of a collapsing sphere, a=2 for early stages and 
a = 1.5 after a si ngular ity form ed at t he centre of the sphere 
(Shu ll977l Larson IT969 1 Penston [19691 Whitworth & Summers 
11985b . Our observed values from the dust (ac=l.9) and C^^O 
(Qfc=L8) PDF slopes are consistent with the spherical free-fall 
scenario, but CS and N 2 H^ have higher values with ac=2.3 and 
ac=2.4, respectively. All high-density clumps traced in CS and 
N 2 H^ emission are associated with pre- and protostellar sources 
(Motte et al. 2007), which indicates that gravitational collapse 
in the star-forming phase has started. The high column densi¬ 
ties seen in dust, CS, and N 2 H^ are thus not a consequence of 
a long column of diffuse emission, but correspond to spatially 
concentrated high volume densities. 

On the other hand, the purely spherical free-fall picture 
seems unlikely to apply to most of the gas which defines the N- 
PDF power laws, because most of the gas we observe in Cygnus 
X is filamentary. This is clearly illustrated in the maps shown 
in Figs. 1, 2, and 7, which display the filamentary network in 
the CO lines and, with only a small fraction of gas in cores, that 
could be considered roughly spherical (see CS and N 2 H^ maps 
in Figs. 2 and 7 and submm continuum observations in Motte et 
al. 2007^ Filaments can be in global free-fall collapse, as sug¬ 
gested by molecular line observations (e.g. Peretto et al. 120131 
Schneider et al. 12010112015bb but the gas we see in filaments 
(and cores) can account for power-law N-PDFs if it has the right 
power-law radial-density structure, even in hydrostatic balance 


density regime. They are thus more precise for separating individual 
clumps. 

^ In Aquila (Konyves et al. 12015b the power-law tail contains more 
than -50% of the mass in filaments but only -15% in cores. 


with no net inward motion. The observed N-PDF slopes (s=- 
1.4 to -2.6, leading to af between 1.4 and 1.7) are consistent 
with s elf-gravitatin g but n on-collapsing filament models (Myers 
120151 Toci & Galli 120151) . In this case, filamentary gas that has 
a sufficient mass per length ratio is centrally concentrated and 
self-gravitating, but does not need to be globally collapsing. 
However, to confine the filaments/clumps, an external pressure 
Pext ^ {^) Thin of the order of pext ~ cm“^ x 10 K is re¬ 
quired (where (n) is the average density in the filaments/clumps 
and Tkin the gas temperature). However, this pressure cannot be 
isotropic, provided by the lower density envelope that is well 
traced in at a density of -10^ cm“^ and a temperature of 
-10-40 K. Moreover, it could be ram pressure of new gas that is 
accreting onto filaments and clumps, a process that is now fre¬ 
quently observe d (e.g . Schneider et al. 120101 Kirk et al. 120131 
Palmeirim et al. 12013b . 

The flatter power-law slopes of the CS and N 2 H^ PDF in 
comparison to the dust could reflect an abundance change be¬ 
cause N 2 H^ is produced wh en CO is depleted at very high col¬ 
umn densities (Bergin et al. 120011 Tafalla et al. 2002, 2006). In 
addition, because CS is not thermalized at the lowest column 
densities, this could lead to a somewhat flatter tail with respect 
to the dust. 

In summary, a power-law distribution in the density struc¬ 
ture is required to have a power-law tail in the PDF and this 
can be achieved by either a hydrostatic configuration, where the 
power law arises from a balance of gravitational forces and pres¬ 
sure gradients, or directly in a dynamically collapsing system. 
Probably both play a role, with filaments possibly balanced or at 
least contraction-slowed by (magnetic) pressure gradients on the 
one hand and very dense collapsing cores on the other. 

6. Summary 

We derived probability distribution functions (N-PDFs) of H 2 
column density for the Cygnus X North region from dust, ^^CO, 
i^CO, 1^0, CS 2^1, and N 2 H+ 1^0. The determination 
of the H 2 column density from the molecular lines is based on 
standard procedures and abundances, assuming a common exci¬ 
tation temperature of 10 K and LTE. Our findings are: 

• The PDFs of dust and CO isotopologues can be described by 
a lognormal distribution for the low column-density range be¬ 
tween Ay-1 to -15, with a common peak around Ay-5. Though 
line-of-sight contamination and variations in abundance and ex¬ 
citation temperatures introduce large uncertainties in the CO- 
PDFs, the overall correspondence between these observational 
PDFs and the ones from recent simulations (Sziics et al. 2015) is 
good. 

• Optical depth effects are the main reason for the different PDF 
shapes of ^^CO, ^^CO, and for higher column densities. 
The PDF from is cut above Ay-30 because the line be¬ 
comes optically thick. The marginally optically thin line 
shows excess in its PDF for Ay=15-40, while the PDF for the 
optically thin line displays a power-law tail in the same 
Ay=15-40 range. 

• Neither selective photodissociation nor fractionation seem to 
play a significant role for the CO abundances for low column 
densities. 

• Depletion of all CO isotopologues probably sets in for Ay 
greater than around 40 where all PDFs fall off. 

• The PDFs of CS and N 2 H^ only consist of a power-law tail 
that covers a high column density range (Ay- 15 to a few 100 
Ay). Using CS as an example, we discuss the potential influence 
of abundance and excitation temperature variations, and show 
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that these could shift the entire PDF by more than one magni¬ 
tude in column density. To better constrain this uncertainty, we 
‘calibrate’ the molecular line PDF using the dust-emission map 
(i.e. shifting the distribution at a given excitation temperature so 
that it corresponds to the dust PDF). For CS, we obtain an abun¬ 
dance of [CS]/[H 2 ] = 10“^ in this way. 

• We find that the N 2 H^ and CS lines are mainly thermally 
excited. N 2 H^ is well mixed with the dust and traces spatially 
the same dense gas clumps of a typical size of ~0.6-0.7 pc that 
were shown to contain pre- and protostellar sources (Motte et al. 
2007). 

• The slopes of the high column-density power-law tail end of 
the N-PDFs from dust, CS, and N 2 H^ are s=-23, -1.6, and -1.4, 
respectively. These values are consistent with a gas distribution 
that is dominated by gravity, i.e. free-falling gas in cores and fil¬ 
aments, although a hydrostatic configuration with ram pressure 
by gas accretion can also take place. The power law then arises 
from a balance of gravitational forces and pressure gradients. 

Summarising the observational results, we find that the N- 
PDFs obtained from molecular lines are not well confined be¬ 
cause they depend strongly on excitation temperature and abun¬ 
dance, and various combinations of these can lead to the same 
PDF. ‘Calibrating’ a molecular line PDF with dust is an appeal¬ 
ing approach but it should be kept in mind that dust PDFs also 
suffer from uncertainties. For example, the specific dust opac¬ 
ity is not constant (a value of 13=2 was chosen for this paper) 
and line-of-sight contamination can lead to an overestimation of 
the column density. However, this study nevertheless shows that 
the dust provides the highest dynamic range in tracing the low 
to high column-density regime, compared to the molecular line 
tracers (see also Goodman et al. 12009 1 and Burkhart et al. 120131) . 
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Appendix A: Determination of H 2 coiumn density 
from dust 

We assume optically thin dust emission at one dust temperature 
Td for each pixel in the map and fitted a greybody function of 
the form ly = By{Td)Kl. to the absolutely calibrated datapoints 
at 160 yum, 250 yum, 350 yum, and 500 yum. The specific dust 
opacity per unit mass (dust+gas) is approximated by the power 
law Ky = 0.1 (y/lOOOGHz)^ cm^g with 13=2, and the dust tem¬ 
perature and surface density distribution E are left as free pa¬ 
rameters. The H 2 column density Ndust is then calculated from 
= fiH niH Ndust ^ adopting a mean molecular weight per hydro¬ 
gen molecule yUj^=2.8. For more details, see Hill et al. (2011). We 
approximate the final uncertainties in the dust column-density 
maps to be around -30-50%, mainly as a result of the uncer¬ 
tainty in the assumed form of the opacity law and possible tem¬ 
perature gradients along the line of sight. 

Appendix B: Determination of H 2 coiumn density 
from moiecuiar iines 


Tex from 1-0 cl= -4 km/s 



Fig.B.l. Map of excitation temperature of the main velocity 
component in Cygnus X North obtained from a Gaussian fit to 
the 1^0 line. 


B.1. Determination of and constraints on the excitation 
temperature 

The temperature stucture in the Cygnus X North region is com¬ 
plex. Large-scale UV illumination from the central Cyg OB2 
cluster heats the diffuse gas and PDR surfaces. Internal (proto)- 
stars and clusters such as DR21 show up as local hotspots. In 
parallel, the Cygnus X North clouds also contain a large reser¬ 
voir of cold and dense gas, seen in submm continuum (Motte et 
al.[2002l). 


The observations of several molecular lines enable us 
to trace sensitively this temperature structure. If the LTE 
assumption is valid and the density is larger than the critical 
density of the transition, the kinetic temperature of the gas (Ty^/„) 
corresponds to the excitation temperature Tex for the molecular 
lines (which is assumed to be equal for all species). If gas and 
dust are well mixed, this temperature should also correspond 
to the dust temperature Tdust that was derived from Herschel. 
However, LTE conditions are not present and the density varies 
from low values in the diffuse gas phase, mainly seen in ^^CO, 
to high values in dense clumps in the molecular clouds, traced 
by dust and the CS and N 2 H^ lines. In the following, we 
determine the excitation temperature in different ways, using 
the observational data sets. 

Tex from 1^0 

Erom 1^0, we calculate the excitation temperature assum¬ 
ing that this line is optically thick so that 

Tex( CO) - + 5.53/(T„i,(i2CO) + 0.818))‘ 

The molecular clouds in the Cygnus region consist of several 
velocity components (Schneider et al. 12006b . We performed a 
Gaussian line fit to the four main l ines (see explanations for the 
column density). Eigure |BT] shows the resulting excitation 
temperature for the most prominent component at velocities 
around -4 km s“^ It becomes obvious that the temperature 
varies significantly between ~5 and -40 K. 

Tex from N 2 H+ 1^0 

We calculate the excitation temperature from N 2 H^ by fitting 
its hyperfine structure with the known relative intensities of the 
seven components. The fit delivers simultaneously the opacity 
T and Tex and shows no anomaly in the relative intensities. 
The average value of Tex is 7 K (in a range of 5 to 20 K) and 
the opacity is generally below one. However, we emphasise 
that only the product rxTex is constrained by the fit so that the 
excitation temperature is not well constrained. 

Tex from dust 

The right-hand panel of Eig. 6 shows the dust temperature we 
derive from the SED fit to the Herschel fluxes. The average 
value for the dust temperature across the map is 15 K, with 
a variation between 10 and 25 K (see also Hennemann et al. 
2012). Because the dust temperature is an average along the 
line of sight, it is not well determined in regions with strong 
temperature gradients. However, for cold, dense clumps and 
cores that are associated with CS and N 2 H^ emission, we expect 
a much smaller variation in temperature. 

In summary, considering the excitation temperatures deter¬ 
mined from the molecular lines and the dust, we use 10 K in the 
following as the best compromise for the determination of col¬ 
umn densities. However, we also discuss other methods (using 
Tex from ^^CO) to obtain the column density, and study 
the influence of different Tex for the N 2 H^ column density. 


B.2. Column densities from linear molecules 


The beam-averaged total column density N of any optically thin 
molecule can be determined from the observed line-integrated 
main-beam brightness temperature with 


N[cm-2] = f(Tex) I T^i,mdv[kms 


f- 
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(B.2) 


and for linear molecules 


f(Tex) = 


3hZ exp(Eup/kTe;,) 

[1 - cxp(-hv/kTeMJ(Te,) - J(Tbg)) 

(B.3) 


and 

J(Tex) = jz ..jrj. X 77 (B.4) 

k(Qxp(hy/(kTex) - 1 ) 

and J(Tbg) = J(2JK), and in which h and k denote the Planck 
and the Boltzman constants, respectively, E^p is the energy of 
the upper level, v is the frequency [GHz], yu is the dipole moment 
[Debye], Ju is the upper value of the rotational quantum number 
and f Tynb dy is the velocity integrated line intensity on a main 
beam temperature scale. 

The first two terms of the rotational partition function for a 
diatomic linear molecule which is accurate to 1%, compared to 
the full term for temperatures above 2 K (Mangum & Shirley 
2015), are given by 

z = ^ + 1/3, (B.5) 

hB 

with the rotational constant B expressed to first order as y=2BJj^. 

To determine all H 2 column densities, we apply a correction 
to the hydrogen mass of a factor of 1.36 to account for helium 
and other heavy elements. 


'^CO 1-0 opacity cl= -4 km/s from LTE (T3,,(12C0)) 



2000 1000 0 -1000 -2000 
offset ["] 

Fig. B.3. Map of the opacity derived for the velocity com¬ 
ponent -4 km s“^ White pixels indicate positions that were not 
well-fitted (e.g. FWHM too large, intensities too low) and thus 
attributed a blanking value. 


H 2 column density from 1-0 

The values for the 1^0 transition are y= 110.201 GHz, 
hy/k=5.29 K, yu=0.112 Debye, 4=1, and /(2.7i^)=0.868. The 
temperature-dependent factor f(Tejc) is 3.22, 2.00,1.62,1.58, and 
1.70x10^^ K“^ (km/s)“^ cm“^ for 5, 7, 10, 15, and 20 K, respec¬ 
tively. As discussed above, we use a constant value of 10 K for 
Tex- A change in the excitation temperature between 5 and 15 


K implies only a 10% variation in the column-density value. We 
then use an abundance [^^C]/[^^C]=70, determined from Wilson 
& Rood (11994b . This value corresponds well to the one derived 
for Orion A (Danger & Penzias [T9^ . For the [H^j^^CO] abun¬ 
dance, we use a value of 1.1 10"^ (Pineda et al. 120101 Fontani et 
al. [20T2t . 

Because it is known that the 1^0 line is not always 
optically thin, an opacity correction ti 3/(1 - exp (r^)) for mod¬ 
erate T can be applied for the column density (Frerking et al. 
11982b . To verify if this kind of correction is required, we deter¬ 
mine the opacity by using the excitation temperature Tex 
that was obtained from and the observed main-beam 
brightness temperature T^bC^CO) with 

t ,3 = - ln(l - T„,(‘3CO)/(5.29/(exp(5.29/T,J - 1) - 0.868)). 

(B.6) 


Because the molecular clouds in the Cygnus region consist of 
several velocity features (Schneider et al. 12006b . we perform a 
Gaussian line fit to four components for and (see 
B.l). As shown in Fig. |B.2[ the dominating emission comes from 
the component at -4 km s“^ , which is associated with the molec¬ 
ular clouds of DR21, DR22, and DR23. The component at 9 km 
s“^ becomes important for the W75 region and southwest of the 
DR21 ridge. Emission at 3 and 15 km s“^ is more diffuse and 
on a low intensity level. We calculate the opacity for all 
velocity components and find that the line is optically thin ev¬ 
erywhere, even for the -4 km s“^ component that is shown in 
Fig. |B.3[ We thus did not include an opacity correction for the 
determination of the column density. 

The N-PDFs that were determined from the column 
density, using a constant excitation temperature of 10 K or a 
variable pixel-to-pixel temperature Tex from ^^CO, do not show 
a significant difference. 


H 2 column density from 1-0 

The values for the 1^0 transition are y= 109.782 GHz, 
hy/k=5.27 K, yu=0.1098 Debye, 4=1, and /(2.7i^)=0.872 K. 
The temperature-dependent factor f(Tex) is 3.34, 2.08, 1.69, 
1.65, and 1.78x10^^ (km/s)-^ cm'^ for 5, 7, 10, 15, and 

20 K, respectively. As discussed in Sec. B.l, we adopt a value 
of 10 K for Tex- To determine the H 2 column density, we then 
use a ratio [H 2 ]/[^^CO] of 1.1 10"^ (Pineda et al. 120101 Fontani 
et al. 12012b and an abundance [^^0]/[^^0] = 531, which was de¬ 
termined from Wilson & Rood (11994b . using a distance of 1.4 
kpc. 

We determine the opacity in the same way for the four 
line components as we did for (see above), using the ex¬ 
citation temperature Tex obtained from and the observed 
main-beam brightness temperature T^/,(C^^O) with 

ti8 = - ln(l - T„i,(C‘«0)/(5.27/(exp(5.27/T„) - 1) - 0.87)). 

(B.7) 

Maps of the resulting opacity show that the line is opti¬ 
cally thin for all velocity components and reaches its maximum 
value of ti8 ~0.4 only in the DR21 ridge for the -4 km s“^ 
component. 


H 2 column density from 1^0 

The 1^0 line is optically thick so that it should a priori 
not be a good tracer for the total column density. However, 
empirical studies (Strong et al. 119881 Bolatto et al. l2Q13b show 
that the line still can trace the total column density/mass of 
a molecular cloud reasonably well with a conversion factor 
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Fig.B.2. Maps of the Gaussian fits to the main-beam brightness temperature of 1^0 for the four main velocity ranges in 
Cygnus X North. The most prominent emission comes from the velocity component at -4 km s“^ 


of 2-2.3x10^^ from line integrated CO main-beam brightness 
temperature into H 2 column density. Here, we adopt a value of 
2x10^^ K-i (km/s)-i cm-^. 

H 2 column density from CS 2^1 

The values for the CS 2^1 transition are y=98.0 GHz, 
hy/k=4.70 K, yu=1.95 Debye, 4=2, and /(2.7X)=1.00. The 
factor f(Te^) is 16.1, 9.48, 7.35, 7.00, and 7 . 41 x 10^2 
(km/s)“^ cm“^ for 5, 7, 10, 15, and 20 K. We assume optically 
thin CS emission and an abundance [CS]/[H 2 ] of 4x10“^^ 
(Fechtenbaum et al., in prep.). 

H 2 column density from N 2 H^ 1^0 

The values for the N 2 H^ 1^0 transition are y=93.17 GHz, 


hy/k=4.47 K, yu=3.40 Debye, 4=1, and /(2.7X)=1.06. The fac¬ 
tor f(Tex) is 3.62, 2.44, 2.10, 2.16, and 2.39x10^2 K-^ (km/s)-^ 
cm“^ for 5, 7, 10, 15, and 20 K. The total beam averaged 
N 2 H^ column density is calculated using the line-integrated 
emission from the F=(2,l) J=(1,0) component at a frequency 
of 93.1737767 GHz and a HFS intensity ratio of 7/27 for that 
component. We assume an abundance [N 2 H^]/[H 2 ] of 5x10“^^. 

We emphasise that all H 2 column-density determinations 
have a large uncertainty, mainly arising from the abundance. The 
[C0]/[H2] abundance for the CO isotopolgues is probably bet¬ 
ter constrained because it has been well studied in the literature. 
However, the variation in excitation temperature along the line of 
sight is more significant here because the CO lines trace a larger 
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range in densities and temperatures than the high-density tracers 
CS and N 2 H^. For the latter, the constant excitation tempera¬ 
ture assumption is most likely more valid, but the abundances of 
CS and N 2 H^ are less well constrained, in particular for massive 
clumps and cores. In summary, we consider the H 2 column den¬ 
sities from CO to be correct within a factor of two and the ones 
from CS and N 2 H^ within a factor of a few. 


Appendix C: Individual molecular line PDFs 

In the main body of the paper, we used a normalization for 
the PDF with respect to the size, i.e. pixel number of the dust 
map. We here show the PDFs obtained with a normalization 

X +00 ^+00 

^ Prjdr] = PncIN = 1. As such, the PDFs are consistent 
with what is shown in the literature (see Schneider et al. 2015a 
and references therein) and the number of pixels defining the 
PDF can be extracted. Figure C.l shows the resulting PDFs. 

Appendix D: Relation between PDF slope and radial 
density distribution for different geometries 

We define the PDF as the surface fraction with gas of column 
density N and the PDF of the natural log of the column density 
density with p = \n{N/{N)), the relation between the two is 
Prj = Np^. The definition of p^ gives Pn ^ dA/dN. 

Spherical geometry 

In the case of spherical geometry, the density radial profile is 
p(r) oc r“^, hence N oc pr oc (dN oc r~^). Hence, with the 

area A oc (dA oc r) 

Pn oc dAjdN oc rlr~^ oc 

p^ = Npn (D.l) 

Cylindrical geometry 

In the case of cylindrical geometry, the density radial profile 
from the axis of the cylinder is p(rc) oc r“^, hence N oc pvc oc 
However in this case A oc ^ x rc oc rc, hence 

Pn oc dAjdN oc llr~^ oc 

p„ = NPN oc (D.2) 

Note that similar calculations can be found in Feder rath & 
Kies sen (12013b . Fischera et al. (12014b . and Myers (12015b . 
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Fig. C.l. PDFs derived from the molecular line maps. The left y-axis gives the normalized probability p(rj), the right y-axis the 
number of pixels per log bin. The upper x-axis is the visual extinction and the lower x-axis the logarithm of the normalized column 
density. For the CO lines, we fitted a lognormal to the low column-density range, indicated by a green curve. The red line shows the 
power-law fit to the high column-density end for C^^O, CS, and N 2 H^. The respective slopes 5' (and errors) are given in the panels. 
For the CO lines, the peak and width of the PDF is indicated with Ay peak and cr. 


13 













